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Короткий зміст 


Ө Організація курсу 


Ө Основи інвестиційного ринку 
© Акції 
ә Біржові товари 
ә Облігації 
ә Ринок Рогех 


Ө вартість грошей в часі 
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Організація курсу 


е Лекційний курс 

е Захист практичних завдань - 80 балів 

ә Теоретичний тест за матеріалами лекцій - 20 балів 
ә Залік 
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Короткий зміст курсу 


Основи інвестиційного ринку. Вартість грошей у часі. Акції 
та облігації. Біржові товари 

Стохастичні процеси. Вінерівські процеси. Броунівський 
рух. Геометричний броунівський рух. Процес Ito. Лема Іто. 
Метод Ейлера-Маруями для розв'язування стохастичних 
диференціальних рівнянь Іто 

Методи Мільштейна та Рунге-Кутти для розв'язування 
стохастичних диференціальних рівнянь [то 

Опціони. Call/Put опціони. Європейський, американський 
та екзотичні опціони 

Модель Блека-Шоулза (Black-Scholes model). Виведення 
рівнянь. Приведення до рівняння дифузії. Аналітичний 
розв'язок для європейського опціону 

Модель Блека-Шоулза (Black-Scholes model). The greeks 
(delta, gamma, vega, theta, rho). Метод Монте-Карло для 
оцінки європейського опціону 
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Короткий зміст курсу 


Метод скінченних різниць для моделі Блека- Шоулза 
Метод скінченних елементів для моделі Блека-Шоулза 


Біноміальні дерева для оцінки вартості опціону 


Моделювання відсоткових ставок. Процес 
Орнштейна-Уленбека. Модель Васічека 

Облігації з нульовим купоном (дисконтні облігації). 
Застосування моделі Васічека для оцінки їх ціни 
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Основи інвестиційного ринку 


ә Акції 


Облігації 


Біржові товари 


Ринок Forex 
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Акції (Stocks) 


ә Акції відображають володіння частиною компанії 


ә Акції надають власнику право на відповідну пропорцію 
майна та доходу компанії 


ə Share (доля акції) - найменша частина акції компанії, що 
може бути придбана 


® Компанії можуть випускати акції, щоб інвестори 
профінансували втілення в життя певної ідеї, яка згодом 
принесе доходи 


ә Акції вважаються ризикованою інвестицією 


ө Власники акцій можуть отримувати дивіденди 
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Фондові біржі та брокерські фірми(посередники) 


ə NASDAQ - National Association of Securities Dealers 
Automated Quotation 

ә NYSE - New York Stock Exchange 

ə NYSE American, колишня American Stock Exchange (AMEX) 


ЕЕ 


FinanceHome Watchlists | МуРогіоїо  Cryptocurrencies Screeners Markets News PersonalFinance Videos | У 


sap 500 Dow30 Nasdaq Russell 2000 Crude ой 
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ssam EAA otoa МО anog S aap ЧОЕ | озо(дзы 


open паз заз 
ва атаах pe satio tiru) asas 
as a7oa5x900 єртї) вз 
бауздаме 17131-17974 ып оазе 202. 


52 мев валд 12904-18284 092 (оа) 
ká avia niat 


Volume 70346295 хе Аца 05,2022 


Avg. Volume 76263275 18376 
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Біржові товари (Commodities) 


ә Commodities - сировинні товари, такі як нафта, природний 
газ, золото 


ә Робити інвестицію в commodities не так просто напряму 


ә Тому виникає поняття ф'ючерсів (future contracts) - 
контракт на купівлю певної сировини в певний час в 
майбутньому за ціною, зафіксованою на даний момент 


ә Також ф'ючерси дозволяють компаніям, що залежать від 
поставок сировини, згладити вплив коливання цін на ці 
товари на дохідність компанії 

ә Володіння сировинними товарами може допомогти 
згладити вплив інфляції, так як ціни на ці товари зазвичай 
зростають при значному збільшенні рівня інфляції 
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Облігації (Bonds) 


ә Облігація - це боргове зобов'язання, яке означає, що 
інвестор позичає гроші комусь (компанії, уряду) 

ə Гроші даються в позику на певний наперед визначений 
період часу та під певну фіксовану чи плаваючу процентну 
ставку 


® Коли компанія хоче профінансувати свої нові проекти вона 
може випустити облігації замість того, щоб брати кредити 
в банках 


ә Відсоткові ставки за облігаціями є зазвичай дещо 
більшими, ніж відсоткові ставки за депозитами банків 


ә Ризики дефолту компанії чи уряду, що випустили облігації 


Віталій Стельмащук Комп'ютерні методи фінансової математики 12 / 15 


Ринок Forex 


ә Forex (Foreign Exchange) - біржа обміну валют 


® Теж посередниками можуть виступати брокерські фірми 
ө Учасниками ринку FOREX є: 
ә Центральні та комерційні банки; 
ә Компанії, що здійснюють зовнішньоторговельні операції; 
Інвестиційні фонди та організації; 
Валютні біржі; 
Брокерські компанії; 
Дилінгові центри; 
Приватні особи. 
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Вартість грошей в часі 


ә Вартість суми х грошей сьогодні більша за вартість тієї ж 
суми х грошей через рік, бо, наприклад, інвестор може 
купити облігацій на цю суму грошей або ж покласти їх на 
депозит 

ə Майбутня вартість FV (Future Value) - скільки буде 
коштувати товар, що зараз коштує х у майбутньому 
(наприклад, через п років) 

ә Теперішня вартість РМ (Present Value) - яка сума грошей у 
даний момент відповідає сумі х у майбутньому 
(наприклад, через п років) 
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Вартість грошей в часі. Дискретна та неперевна моделі 


ә Дискретна модель. Вважаємо, що ринкова відсоткова 
ставка г - стала величина. Тоді 


РУ = хх (1+ г)", 
п a (1) 
(1+) › 
де п - кількість років. 
ә Неперервна модель. Ринкова відсоткова ставка г - стала 
величина. Ми хочемо знати, яка вартість грошей x(t) в 
кожен момент часу. Тоді 


x(t + dt) — x(t) = Фа = rx(t)dt, 


& = i(t), (2) 
x(t) = х(0)е", 


тобто 
FV(t) = PV x e*, (3) 
PV = FV (t) x e7". 
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Стохастичні процеси. Вінерівські процеси. 
Броунівський рух. Лема |то. Геометричний 


броунівський рух 
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Стохастичні процеси. Базові поняття 


ə Стохастичний (випадковий процес) - сімейство 
випадкових змінних Х;, що визначене відносно системи 
параметрів t. Часто t € [0, Т] - це час, а Х, = X(t) - 
випадкова змінна для кожного t € [0, Т]. 

ә Неперевний в часі стохастичний процес {Х;; t є [0, Т]} 
називається процесом Гауса тоді і тільки тоді, якщо для 
будь-якої скінченної множини індексів Ы,..., їх з множини 
10, Т] випадкова змінна Х,, п = (Хн,..., Хь) має 
багатовимірний нормальний розподіл. Іншими словами: 
кожна лінійна комбінація (Хн,..., Ху,) має одновимірний 
нормальний розподіл. 

ә Неперевний в часі стохастичний процес { Хе; t є |0, Т|) 
називається процесом Маркова (чи має властивість 
Маркова), якщо на його наступне значення впливає лише 
теперішнє значення Xe. 
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Броунівський рух. Вінерівські процеси 


ә Броунівський рух названо на честь Роберта Броуна, який 
ще в 1827 році спостерігав рух частинок пилку в воді. 


ə Вперше математично описав такий рух частинок Норберт 
Вінер у 1923-1924 роках. 


ә Броунвський рух називають ще "випадковим 
блуканням" (гапфот walk), а також вінерівським процесом. 


ә Стохастичний процес We, Е Є [0, T] називають 
вінерівським, якщо: 

ә W; - неперервний i И = 0 

ə W; ~ N(0, t) Yt > 0, тобто для кожного t випадкова змінна 
W; розподілена нормально з математичним сподіванням 
Е(И/,) = 0 та дисперсією М(М/,) = Е(И/2) = t. 

ә Приріст dW; = И/ ++ — Wr - нормально розподілений, 
тобто dW; ~ М(0, dt), а отже Е(АМ/,) = 0 та М(АМУ,) = dt. 
З цього також випливає, що приріст dW; не залежить від 
історії процесу до моменту часу t. 
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Вінерівський процес. Симуляція 


ə Розбиваємо часовий інтервал [0, T] на М підінтервалів, 
ДЕЕТ 
ә Починаємо з Wo = 0 
ә Для кожного j = 1,2..., №: 
возі + At; 
AW ~ УДЕМ(0, 1); 
W; = Wi- + AW. 
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Вінерівський процес. Симуляція 


Улепег-ргосе55 


М/іепег-ргосез5 W(t) 


т T T T 
200 400 600 800 1000 
Time(t) 


o4 


Puc.: Приклад вінерівського процесу 
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Броунівська модель ціни фінансових активів 


е Вінерівський процес може набувати від'ємних значень, 
тому очевидно, що він не дуже підходить до моделювання 
цін на активи. 

ə Дохідність активу (наприклад акції) в момент часу t 
визначається за формулою: 


ЕЕ Зена = 5, Е dS: 
- 5, т 5, И (1) 


де 5; - ціна активу в момент часу Ї. 


К, 


® Дохідність моделюється броунівським рухом з дрейфом 

(drift) i параметром волатильності (volatility): 

dS 

Р, = = udt + odW,, (2) 

5 
де udt - визначає тренд дохідності активу, саМ/, - 
відповідає за випадковість. Тут | - параметр дрейфу (в 
процентах), с - параметр волатильності (в процентах). 
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Броунівська модель ціни фінансових активів. 


Геометричний броунівський рух 


ә З рівняння (2) випливає: 
dS: = uStdt + oStdWe, (3) 


ә Стохастичний процес St, що задовольняє стохастичне 
диференціальне рівняння (3) називають геометричним 
броунівським рухом (Geometric Brownian Motion). 

® Відносні, а не абсолютні зміни в ціні активу 
задовольняють броунівський рух з дрейфом. 
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Процес lto. Лема |то 


Процес Іто 


Стохастичний процес Х;, що задовольняє стохастичне 
диференціальне рівняння Іто 


ах, = а(Х;, гаї =r b(Xt, ї)а\уМ,, (4) 


називають процесом Іто (Ito process). 
а(Ху, t) - коефіцієнт дрейфу (drift), 
Ь(Хь, t) - коефіцієнт дифузії (diffusion), 
ИИ, - вінерівський процес. 


ə Аналог правила диференціювання складної детермінованої 
функції: 


4869,9 


ве a 
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Лема |то 


Нехай Х; задовольняє рівняння процесу Іто: 
ах, = а(Х,, t)dt + ЫХ, Вами, (6) 


та &(х, t) - С?1-гладка функція (існують неперевні похідні 


а Peg дє) 
Ох’ дх? › дЕ/` . 
Тоді випадкова змінна У, := 2(Х;, t) задовольняє процес |то з таким 


самим вінерівським процесом ИЖ: 


а ТЕ. дЕ 
ау, = (Zas зо доб Пол ВО (7) 


де всі похідні від g і коефіцієнти а та b насправді залежать від (X+, t). 


Доведення 


Через розклад функції від двох змінних у ряд Тейлора і оцінку 
кожного доданка при dt — 0. Детальніше, наприклад тут: 

Levy G. Computational Finance Using С апа C# / С. Levy. - Academic 
Dra An 
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Геометричний броунівський рух. Аналітичний розв'язок 


ә Геометричний броунівський рух описується стохастичним 
рівнянням: 


45; = паї + о5,а\М,, (8) 
ө Лема Іто для цього рівняння: 


_ (08 дав 108g 222 де 
ас (5, t) = (2215.4 ді T 2252 S: озна, (9) 


ә Розглянемо функцію 2(5., t) = In(S+). 


ө Тоді 
де _ 9(($.)) _ 1 
б во 
В — = 
95,2 — 95, (2) эш: (10) 
дв _ 
ðt 
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Геометричний броунівський рух. Аналітичний розв'язок 


ә Підставимо похідні (10) у рівняння (9): 


1 
де (5., t) == (в- 2") dt + сам, (11) 
тобто g(St, t) = In(S+) задовольняє рівняння броунівського 


руху з дрейфом (и = 102) та волатильністю с. 


ә Проінтегрувавши ліву і праву частини рівності (11), 
отримаємо: 


1 

п.) = | рео | ее Yt € (0, Т]. (12) 
55 2 
ә З рівності (12) отримаємо: 


5, = буе 292) +0: ур є (0, Т]. (13) 
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Метод Ейлера-Маруями для розв'язування 


стохастичних диференціальних рівнянь [то 
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Апроксимація Ейлера-Маруями 


Стохастичне диференціальне рівняння [то 


АХ, = а(Хь, t)dt + Ь(Ху, t)dW:, x Є По, Т], Xo = Хо (1) 
Апроксимація Ейлера-Маруями 

Неперевний в часі стохастичний процес У = {Y (t), к < І < T} 
називається дискретною апроксимашею Ейлера-Маруями, якщо 
для заданої дискретизації to = To < ту < ... < тм = Т часового 
інтервалу [&, Т] вн задовольняє ітераційну схему: 


Ул = Yo + а( а = Та) aF b( Yn, Tn) (ИИ = 29 
Yn = 0,1, ..., М—1 


1Гулг Yo = Еу 
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Алгоритм Ейлера-Маруями 


Апроксимація Ейлера-Маруями 


На практиці часто використовують рівномірну сітку вузлів на 
часовому інтервалі: 
T — tọ 
At = —— 
N 


Введемо позначення: 


AW, = W, 


прі 


— Wn: 


п 


Тоді схема (2) запишеться у вигляді: 


Үд = Yn + a(Yn, та) Де БУ, та) ДАМУ, 
Уп = 0,1,..., М -1 


Алгоритм Ейлера-Маруями 


ә Yo = Хо. Для кожного j = 0,1,2...,М- 1: 
СЕЧА 
AW ~ У ДЕМО, 1); 

Уна = У + а(ү, в) АЕ БУ, ti) AW. 
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Слабкий і сильний розв'язок 


ә Як відомо, рівняння, що описує геометричний броунівський 


рух 
dS: = uStdt + с5,а%, (6) 
має аналітичний розв'язок: 
Se = Soelt 30°) 9 Yt є (0, Т]. (7) 


ә В залежності від конкретної реалізації вінерівського 
процесу W, розв'язок 5; набуватиме різних значень. 
Такий розв'язок 5; для конкретної реалізації вінерівського 
процесу W; називається траєкторією. 

ә Якщо нам задано конкретний вінерівський процес We, то 
тоді розв'язок 5; називається сильним розв'язком. 

ә Якщо ми вільні у виборі вінерівського процесу We, то тоді 
розв'язок 5, називається слабким розв'язком. У випадку 
слабкого розв'язку нас цікавить лише розподіл St, а не 
його індивідуальна траєкторія. 

ә Для інших процесів [то означення аналогічні. 
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Абсолютна похибка 


ə Нехай ў - чисельна апроксимація процесу |то, що 
задовольняє рівняння (1). Тут А := ДЕ. 


ә Взявши ідентичну траєкторію одинакового вінерівського 
процесу процес W; для стохастичного рівняння і для 
чисельної апроксимації, ми можемо порівнювати значення 
Х; та Ү для будь-якого вузла сітки tj. 

е Абсолютна похибка апроксимації сильного розв'язку 
дорівнює |Хт — yaj: Для іншої траєкторії така абсолютна 
похибка буде іншою. Ми зацікавлені в середньому значенні 
абсолютної похибки для всіх траєкторій. 


о Отже, абсолютна похибка апроксимації Ү сильного 
розв'язку Х, рівняння (1) визначається за формулою: 


е(һ) := Е(|Хт - ҮЛ). (8) 
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Оцінка абсолютної похибки на практиці 


ә Проводимо N експериментів, вибираючи іншу конкретну 
траєкторію вінерівського процесу We. 

® Обчислюємо аналітичний і чисельний розв'язки згідно цієї 
траєкторії процесу We. 

ə Порівнюємо розв'язки і усереднюємо: 


N 
р 1 
ё(һ) := N ` Хт, — ү? xl: (9) 
k=1 
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Сильна та слабка збіжність 


ә Сильна збіжність. у збігається сильно до XT з 
порядком збіжності y > 0, якщо 


е(һ) := E(|Xr - УҢ) = O(h). (10) 
ya збігається сильно, якщо 


lim Е(|Ху - ў) = 0. (11) 
һ—›0 


ə Слабка збіжність. ү? збігається слабко до Хт згідно 
функції є з порядком збіжності В > 0, якщо 


Е(2(Хт)) — Е(&(УТ)) = O(h’). (12) 


ҮР збігається слабко з порядком збіжності В > 0, якщо 
(12) виконується для всіх поліномів g. 
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Сильна та слабка збіжність. Метод Ейлера-Маруями 


Відношення між сильною та слабкою збіжністю 


ә Слабка збіжність. З використання полінома g(x) = х 
випливає збіжність математичних сподівань. 


ә Слабка збіжність. З використання полінома g(x) = х2 


випливає збіжність дисперсій. 
ә З сильної збіжності випливає слабка збіжність за 
поліномом g(x) = х, тобто збіжність математичних 
сподівань. 


Метод Ейлера-Маруями і збіжність 


КІоеаеп Р.Е. Numerical Solutions ої Stochastic Differential Equations 
P.E.Kloeden, E. Platen. - Berlin:Springer-Verlag, 1992. - 668p. 
ə Апроксимації методу Ейлера-Маруями сильно збіжні з 
порядком збіжності y = 1. 
ә Апроксимації методу Ейлера-Маруями слабко збіжні з 
порядком збіжності 5 = 1. 
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Методи Мільштейна та Рунге-Кутти для 


розв'язування стохастичних диференціальних 


рівнянь [то 


Стельмащук Віталій Володимирович 


Львівський національний університет імені Івана Франка 
Кафедра інформаційних систем 


27 жовтня 2022 р. 
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Апроксимація Мільштейна 


Стохастичне диференціальне рівняння [то 


АХ, = (Хо, + ЫХ, а, х є |, гк = Хо (1) 


Апроксимація Мільштейна 


Неперевний в часі стохастичний процес У = {Y (t), к <t< T} 
називається дискретною апроксимацією Мільштейна, якщо для 
заданої дискретизації & = то < ту <... < тм = Т часового 
інтервалу [ю, Т] вн задовольняє ітераційну схему: 


ЖЕЛГЕ ПАТО ВУ АИ 
+4b(Yn, Tn)b' (Yn, Tn) (AWM)? — АЕ (2) 
Уп = 0,1,..., М-1 
Тут Yo = Ху, Yn = Y (Tn), AW, = W, 


п+1 
Виводиться з розкладу в ряд Іто- Тейлора. 


— Wn- 


п 
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Алгоритм Мільштейна 


Алгоритм Мільштейна 


о % = Хо. Аг = 10. 
Для кожного j = 0,1,2..., М- 1: 
да = tj + At; 
AW ~ VAtN(0,1); 
Ура = Y + а(ү, &)ДЕ+ bO AW + 
ЗЫ 5)Ы(Ү», 5) (AWP — А4. 


Сильна та слабка збіжність 


ә Метод Мільштейна сильно збіжний з порядком уу = 1. 
ə Метод Мільштейна слабко збіжний з порядком beta = 1. 


ə Доведення: Kloeden Р.Е. Numerical Solutions ої Stochastic 
Differential Equations P.E.Kloeden, E. Platen. - 
Berlin:Springer-Verlag, 1992. - 668p. 

ә Коли b(X:, t) насправді не залежить від Xt, метод 

Мільштейна вироджується в метод Ейлера-Маруями. 


Віталій Стельмащук 3/6 


Схеми з сильною збіжністю вищих порядків 


Схема Платена-Вагнера. Порядок збіжності 1.5 
(4.1) У = У +aA+bAW + Eo (ди)? - А} 
П 1 t 1 2” 2 
+аФА2+ 5 аа +58 a”) А 
+ (w + zee) {AW A – AZ} 
1 и n2 1 2. 
+5 (+) ) {jam д{ ди. 
Неге Ше additional random variable AZ із required to represent the double 


integral 


Тачі раз 
(42) AZ = lag = / / dW,, зэ. 
Рис.: Схема Платена-Вагнера 


Kloeden P.E. Numerical Solutions ої Stochastic Differential Equations 
P.E.Kloeden, E. Platen. - Berlin:Springer-Verlag, 1992. - 668p. 
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Методи Рунге-Кутти 


Недоліком методів побудованих за допомогою розвинень в 
ряд Тейлора є необхідність обчислення похідних а’, Б),.... 


По аналогії з випадком звичайних диференціальних рівнянь 
існують альтернативні методи, які вимагають тільки 
обчислень значень функцій а, б, а не їх похідних різних 
порядків. Такі методи називають методами Рунге-Кутти. 


З методу Мільштейна виведемо альтернативний метод 
Рунге-Кутти з тим же ж порядком сильної збіжності. 


Почнемо з рівності: 


Бу + Ду) - Бу) = Б(у)Ду + О((Ду)?) (3) 
Враховуючи, що Ау = аА? + БДИ, отримаємо: 


Му + Ду) - Бу) = b'(y)(aAt + БАМ/) + О((Ду)?) (4) 
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Метод Рунге-Кутти з сильною збіжністю порядка y = 1 


ә Оскільки Е((ДИ/)?) = At, то рівність (4) можна 
спростити до вигляду: 
Бу + Ду) - Ку) = b'(y)b(y)AW + О(А) (5) 
ә Підставимо AW = у At i отримаємо таку оцінку добутку 


b(y)b'(y) з точністю О(//Дт): 


Бу) (у) = —== (Му + aly)At + bly) АЮ) - Бу). (6) 


1 
УДЕ 
ә В результаті отримуємо таку схему методу Рунге-Кутти: 


ӯ зу; + адёЕ+ БУ А 
уна = yj +aAt+ ИИ + АИ” — At) [0(9) - Б(у)) 


(7) 
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Опціони. Са|/Ри{ опціони. Європейський, 


американський та екзотичні опціони 


Стельмащук Віталій Володимирович 
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і 


Деривативи (Derivatives) 


ə Дериватив - це укладений контракт щодо якогось 
базового активу (underlying asset). Ціна дериватива 
визначається виходячи з ціни базового активу. 


е Приклади базового активу: акції, облігації, біржові товари, 
відсоткові ставки. 
ə Використовуються для: 
ә Обмеження ризику (hedge risk). 


Приклад. Ціни на нафту i відповідно на пальне - дуже волатильні. Компанії, 
чия діяльність надто залежить від пального (наприклад, авіакомпанії) дуже 
ризикують. Тому вони укладають договори з нафтовидобувними компаніями про 
закупівлю нафти/ пального за фіксованою наперед ціною через деякий час у 
майбутньому (через декілька місяців, рік, тощо). Таким чином вони страхують себе 
від потенційних збитків, якщо ціна на нафту сильно зросте. З іншого боку, такий 
договір також вигідний нафтовидобувним компаніям, адже вони страхують себе від 
стрімкого падіння цін на нафту. 

ə Спекуляції (отримати вигоду від того, що ціна на базовий 
актив зміниться в очікуваному інвестором напрямку). 
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Типи деривативів 


Форвардний контракт - це договір укладений між двома 
сторонами про купівлю/продажу базового активу (акції, 
облігації, біржові товари) за наперед домовленою ціною в 
певний момент у майбутньому. Договір не 
стандартизований і ним не можна торгувати на біржі. 
Ф'ючерсний контракт - такий самий зміст договору, як 
при форвардному контракті, але він стандартизований і 
може торгуватися на біржі. Покупець ф'ючерсного 
контракту має право продати такий контракт протягом 
строку його дії іншим особам без погодження умов такого 
продажу з продавцем контракту. 

Своп (Swap) - договір про обмін товарами в майбутньому 
між сторонами. 

Опціон (Option) - контракт, що дає право (не 
обов'язок) купити чи продати базовий актив за заданою 
ціною в майбутньому. Премія опціону -- це сума грошей, 
що сплачується покупцем опціону продавцеві 
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ә Купимо актив, якщо ціна на нього піднялася: 


s(t) 


Puc.: Ціна активу зросла 


ə Не будемо купувати актив, якщо ціна на нього впала: 
s(t) 


| + + >t 
t=0 т 
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Call опціон 


Call опціон - право купити певний актив (акцію) за 
домовленою ціною Е в момент часу Т в майбутньому. 
ә У випадку Call опціону ми спекулюємо на очікуванні, що 
ціна на актив зросте. 


ә Е - Exercise price (strike price) 
ə T - Expiration date (expiry) / Дата виконання 


ə S(t) - Ціна активу (акції) 
ә V(S,t) - Ціна опціону 
ə М(5,Т) = тах($(Т) - Е, 0) - Ціна опціону в дату 


виконання 
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Риї опціон - право продати певний актив (акцію) за 
домовленою ціною Е в момент часу Т в майбутньому. 
ө У випадку Риї опціону ми спекулюємо на очікуванні, що 
ціна на актив впаде. 


ә Е - Exercise price (strike price) 


T - Expiration date (expiry) / Дата виконання 

ə S(t) - Ціна активу (акції) 

V(S,t) - Ціна опціону 

V(S,T) = max(E-S(T), 0) - Ціна опціону в дату виконання 
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Опціони. Короткі та довгі позиції 


Емітент опціону є короткою стороною (англ. short) контракту, 
власник -- довгою стороною (англ. long). Відповідно є чотири 
можливі позиції в опціонах: 


ә купити купівельний опціон (long call) 
ә емітувати купівельний опціон (short call) 
е купити опціон продажу (long put) 


ә емітувати опціон продажу (short put) 
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Long call. Діаграма прибутку 


S(t) = E 


^ 5(0) 


Puc.: Long call. Діаграма прибутку 
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Long put. Діаграма прибутку 


Puc.: Long put. Діаграма прибутку 
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опейський та американський опціони 


Європейський опціон 


Європейський опціон - це опціон, який може бути виконаний 
тільки в дату виконання Т. 


Американський опціон 


Американський опціон - це опціон, який може бути виконаний в 
будь-який момент часу до дати виконання Т. 


Відношення між американськими та європейськими опціонами 
ə Ціна американського Call опціону завжди дорівнює ціні 
європейського Call опціону 
ə Ціна американського Put опціону зазвичай відрізняється 
(більша) за ціну європейського Put опціону 


[| Ruppert D. Statistics апа Finance. Ап Introduction / D. 
Ruppert. - New York: Springer, 2004. - 496p. 
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Put/Call паритет для європейського опціону 


Table 11.2 Portfolios illustrating put—call parity. 


Portfolio A Call option Sr- К 0 
Zero-coupon bond K K 
Total Sr K 

Portfolio C Put Option 0 K – Sr 
Share Sr Sr 
Total Sr K 


Puc.: Портфелі, що ілюструють Put/Call паритет 


[] Hull J. Options, Futures апа Other Derivatives. 11th ed. / J. 
Hull. - Harlow:Pearson, 2021. - 882p. 
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Put/Call паритет 


Європейський опціон 


С+ Ке" = Р + 5) 


Американський опціон 


АРА 


[] Hull J. Options, Futures апа Other Derivatives. 11th ed. / J. 
Hull. - Harlow:Pearson, 2021. - 882p. 
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Екзотичні опціони 


ә Велика різноманітність 


ә Емітент намагається привабити покупця покращеними 
умовами 
ә Нестандартні американські опціони: 
е Виконати опціон можна лише в певні конкретні дати перед 
датою виконання опціону (Бермудський опціон) 
ə Виконати опціон можна лише в певний період життя 
опціону 
ə Execution price (strike price) може змінюватися під час 
життя опціону 
ө Азійські опціони - опціони, в яких виплата залежить від 
середнього арифметичного ціни базового активу протягом 
життя опціону. Тобто, виплата дорівнює 
тах(0, Saverage — К) для call опціону і тах(0, К — Saverage) 
для риї опціону. 
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Модель Блека-Шоулза (Black-Scholes model). 
Виведення рівнянь. Приведення до рівняння 


дифузії. Аналітичний розв'язок для 


європейського опціону 
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Модель Блека-Шоулза 


Сформульована Блеком і Шоулзом у 1973р. 


1997р - Нобелівська премія за це дослідження 


Дозволяє обчислити ціну опціону У(5, Е) у будь-який 
момент часу t < T 


ə Показує, що комбінування ризикованими активами 
дозволяє уникнути ризику взагалі 
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Модель Блека-Шоулза. Фундаментальні припущення 


Припущення про активи 


® Зміна ціни акції в часі описується геометричним 
броунівським рухом 


ə Ринкова відсоткова ставка за безризиковим активом г - 
стала величина 


ə Дивіденди за акцією не виплачуються 


Припущення про ринок 


ə Немає можливості для арбітражу (безризикового 
спекулятивного доходу) 


® Можливість купувати і продавати будь-яку кількість акцій 
(навіть дробову) 


ә Всі транзакції з активами - безкоштовні (немає комісій) 


ә Торгівля активами - неперервна 
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Модель Блека-Шоулза. Виведення рівняння 


Специфічний портфель 


ə Геометричний броунівський рух для ціни акції: 
45; == и5;аї ав с5.аМИ, (1) 
ә Лема |то для ціни опціону М(5$+, t): 


ду И о дм 
ам (5, Е) = (из: + ді ар 2с? S ) о 
(2) 


ә Ідея в тому, щоб скомбінувати в один портфель опціони на 
акцію і самі акції таким чином, щоб вилучити взагалі 
доданки, що залежать від випадкового вінерівського 
процесу W: 


ə Тому розглянемо портфель, де є —1 опціон (коротка 
позиція) та А акцій (довга позиція). 
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Модель Блека-Шоулза. Виведення рівняння 


ә Вартість портфоліо: 
т(ї) = А5, – (5, #) (3) 
ә Зміна вартості портфоліо dr = r(t + dt) - п(П): 
ат = AdS: — dV (S, t) (4) 
ə Зміна вартості портфоліо після підстановки: 
dr = А (uSidt + с5,аМИ,) 


5 
(205,5 +1025) dt + 2075,40) (5) 
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Модель Блека-Шоулза. Виведення рівняння 


ө Зміна вартості портфоліо після перетворень: 


ат = [(А— BE) uS: - (GE + ау S) dt 


(6) 
-(а - 93) osam, 
ә Якщо виберемо А = о то нам вдасться позбутися 
доданку з вінерівським процесом: 
ду 1ё\ 
dr = | 2S? ) dt 7 
т (9 55-20 2) (7) 
ә З іншого боку, зміна в вартості портфелю не може бути 
більшою, ніж приріст безризикого активу, тому: 
dr = rrdt (8) 


е Прирівнюючи (7) та (8), отримаємо: 


ƏV 18V 3.3 Əv 
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Модель Блека-Шоулза. Виведення рівняння 


ə Скоротивши на dt та позначивши 5 = 5, отримаємо: 


ду ыг дм М 
ато Кы аа. (10) 


ә Рівняння (10) відносно ціни опціона V називається 
рівнянням Блека- Шоулза. 


ә Рівняння (10) - параболічне та зводиться до рівняння 
дифузії (теплопровідності). 

® Для знаходження конкретного розв'язку рівняння (10) 
повинно бути доповнене відповідними кінцевою (в точці 
t = Т) та крайовими умовами. 


ә Для європейського Call опціону ці умови матимуть вигляд: 


V(0,t)=0 Yt <= T 
V(S, t) = S — Ке T-t) Yt <= T, S > оо (11) 
V(S, Т) = тах(5 — K,0) YS 
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ння рівняння Блека-Шоулза до рівняння 


Віталій Стельмащук 


ення рівняння Блека-Шоулза до рівняння 


опровідності 


Рівняння Блека-Шоулза після підстановки (12)-(17) у (10) 


vr = Ухх + (k - Тм, — КУ, 


Ще одна заміна змінних 


р 


v(x, T) = КО) 


vx = 698" (au(x, т) + ux) 
Е (œ? u(x, T) +2aux + uxx) 
vr = e®*+tbT(Bu(x, т) + ur) 
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ення рівняння Блека-Шоулза до рівняння 


лопровідності 


Підстановка (20)-(21) у (18) 


шь = Б о k = Pu 4 0 4 k а. 


Вибіраї В 
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Зведення рівняння Блека-Шоулза до рівняння 
теплопровідності 


2 
иг = Uxx, ХЕ ( - 00,00), т Є (0, и 
и(х, 0) = (е'*+1)х/2 — „(К—1)х/2\+ при x Є (—оо, 00); 


То? 
и(х,т) — 0 при x > —œ для т Є (0, = (25) 
м (e нина £ е 66-1072) 


То? 
при х — оо для т Є (0,2). 
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яння Блека-Шоулза. Аналітичний розв'язок 


Метод розв'язування 


Застосовується перетворення Фур'є до рівняння (25). 
Повний процес розв'язування можна переглянути в 
Buchanan К. Ап Undergraduate Introduction to Financial 
Mathematics. 3rd ed. / R. Buchanan. - Singapore: World Scientific 
Publishing Co. Pte. Ltd., 2012. - 486p. 


Розв'язок. Формули Блека-Шоулза 


С(5, t) = 5М(ф) – Ке "(7 9 Мф) (26) 
P(S, t) = - 5М(-аї ) + Ке TÌ N(—a2) 


Використані позначення 
ОСЕ АКИО 
ovyT -t 

d =dı—ovT -t 
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а 


Формули Блека-Шоулза. Кумулятивна функція 


стандартного нормального розподілу 


Кумулятивна функція стандартного нормального розподілу 


Figure 15.3 Shaded area represents N(x). 


Puc.: Кумулятивна функція стандартного нормального розподілу 


П Hull J. Options, Futures апа Other Derivatives. 11th еа. / J. Hull. – Найом:Реагзоп- 2021. = 882p- 
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Модель Блека-Шоулза (Black-Scholes model). 
The greeks (delta, gamma, vega, theta, rho). 


Метод Монте-Карло для оцінки європейського 
опціону 
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Нехай ціна на акцію рівна 1005, а ціна на відповідний 
са!|-опщон - 10$ 
oV 


Нехай А = — = 0.6. 
ехай 25 


Інвестор продав 2000 саї-опціонів на акцію (short position). 


Тоді щоб захистити себе від ризику втрати, йому необхідно 
купити 0.6 ж 2000 = 1200 акцій (long position). 

Якщо ціна на акцію виросте на 15, то ціна на відповідний 
саії-опціон виросте на 0.6 * 1 = 0.65. Тоді з акцій інвестор 
отримає 12005 прибутку, а з опціонів - 0.692000 = 12005 
втрат. Тому його позиція не зміниться. 

Якщо ж ціна на акцію впаде на 15, то ціна на відповідний 
саії-опціон впаде на 0.6 * 1 = 0.65. Тоді з акцій інвестор 
отримає 1200$ втрат, а з опціонів - 0.672000 = 12005 
прибутку. Тому його позиція знову ж не зміниться. 
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oV З : А : | | 
ө: A= 95 ` чутливість ціни опціону до ціни на відповідну 


акцію. А портфелю - це сума А кожної позиції в портфелі. 
Використовується в A-hedging. 
ФМ . : | В | 

е: Г = 952 7 чутливість А до ціни на відповідну акцію. Г - це 
міра того, як часто портфель має бути re-hedged, щоб він 
залишався в Д-нейтральній позиції. Д-нейтральна позиція 
- це позиція, при якій вартість портфелю залишається 
незмінною при незначних змінах в ціні відповідної акції. 

ә О = - - швидкість зміни вартості опціону в часі. 


дї 


ә Vega = ВЯ швидкість зміни вартості опціону по 
С 


відношенню до волатильності с. 
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The greeks. Рівняння Блека-Шоулза 


Рівняння Блека- ДЕЕ, 
ди 1 2 
ОСОМ 2 os? 


oV 


38 (1) 


Рівняння Блека-Шоулза з The greeks 


1 
Ө + 50° 5°Г-+ /5А – ту = 0 
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Метод Монте-Карло 


ә Не існує чіткого означення, система методів 
® Дозволяє розв'язувати детерміністичні проблеми за 
допомогою ймовірностей 


ə Вперше був використаний в 1940-их роках в Manhattan 
Project 


ə Оцінити площу круга, вписаного в квадрат зі стороною а. 


® Для цього будемо випадковим чином генерувати велику 
кількість точок, що попадають в цей квадрат і будемо 
підраховувати скільки з них попадають і в круг 


ә Тоді площу круга можна оцінити таким чином: 


с г points within the circle (3) 
шр = © - 
д total number ої points 
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Метод Монте-Карло для оцінки вартості європейського 


опціону 


ө Зміна ціни на акцію описується геометричним 
броунівським рухом: 


45; = 5: + с5,а%, (4) 
що має аналітичний розв'язок: 
Se = Ѕуе(и— 29) ур є (0, Т]. (5) 


ә Для великої кількості ітерацій / (не менше 10000) рахуємо 
S(T) за формулою (5), використовуючи щораз іншу 
траєкторію процесу Вінера. 

ә Для кожної ітерації / рахуємо відповідну ціну 
європейського опціону. 

Саї-опціон: 

С-(Т) = тах(5(Т) - Е, 0) (6) 
Риф-опціон: 

РКТ) = тах(Е - 5(Т), 0) (7) 
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Метод Монте-Карло для оцінки вартості європейського 


опціону (продовження) 


ә Рахуємо середнє арифметичне значення ціни опціону для 
всіх ітерацій: 


Саї-опціон: 
N 
>. G(T) 
Caverage( Т) == N (8) 
Риф-опціон: 
N 
>. РКТ) 
Paverage( Т) == N (9) 


® Для оцінки вартості опціону в початковий момент часу 
обчислюємо Present Value від отриманого середнього 
значення. 


С(0) = ааа | Т), 


Р(0) - Т ране T (10) 
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Біноміальні методи для оцінки вартості опціону 
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Біноміальний метод 


ə На практиці часто є зацікавленість тільки в значенні 
вартості опціону в даний момент часу У( 90, 0), коли 
відома ціна самої акції 50. Тому немає змісту обраховувати 
значення У(5, t) для всієї області, щоб тільки виділити з 
отриманих результатів V (So, 0). 


ә Відносно мале завдання з обрахування У(5о, 0) може бути 
виконано з допомогою біноміального методу, який 
базується на дерево-подібній сітці вузлів, в кожному з яких 
застосовуються відповідні бінарні правила. Сітка не є 
визначеною наперед, а будується самим біноміальним 
методом. 
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Біноміальний метод. Дискретна модель 


Поділимо часовий інтервал [0, Т] на М рівних підінтервалів і 
введемо позначення: 


ДЕ:= т: 
t; := iAt, i = 0,..., М; (1) 
S = S(t). 


Припущення біноміального методу 


© Ціна 5 під час кожного періоду часу At може змінитися 
тільки або вгору на Su ("up"), або вниз на Sd ("down"), де 
0<а<и 

Ө Ймовірність руху ціни вгору Р(ир) = ві 

© Математичне сподівання та дисперсія 5 відповідають їхнім 


неперервним відповідникам, що обчислені для відсоткової 
ставки безризикового активу г. 
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Біноміальний метод. Дискретна модель 


Рис.: Біноміальний метод. Дискретна модель 


[] Seydel R.U. Tools for Computational Finance. 6th ed. / R.U. 
Seydel . - London: Springer-Verlag, 2017. - 508p. 
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Біноміальний метод. Виведення рівнянь 


® Математичне сподівання для дискретного випадку: 
E(Si+1) = pS;u + (1 - р)5;а (2) 
ә Математичне сподівання для неперервної моделі: 
Е($ 1) = бе" А! (3) 
е Прирівнюємо (2) та (3): 
e"^t = pu + (1 — p)d. (4) 


ө 3 останньої рівності 


р= не (5) 
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Біноміальний метод. Виведення рівнянь 


ә Дисперсія для неперервної моделі: 


ЕВ) = 5 е (6) 
Уаг(5) = Е(57) - (Е(5))" (7) 

З рівнянь (3) та (6) випливає, що дисперсія дорівнює: 
Маг(Зіча) = SPEA Це” А" — 1). (8) 


ә Дисперсія для дискретної моделі: 
а р= Е(5%1) = (Е(5-1))* п № 
= р(5ш)? + (1 - (5) - 52 (ри + (1 - руд) 
е Прирівнявши (8) та (9), отримаємо: 
ТА — ри? + (1 — р)а. (10) 


ә Маємо 2 рівняння (4) та (10) для визначення 3 невідомих 
и, д,р. Щоб остаточно визначити невідомі, виберемо собі 
третє рівняння: 


иа = у (11) 


Біноміальний метод. Модель Кокса-Росса-Рубінштейна 


о У рівнянні (11) вибирається 
тет (12) 
ә Дерево з вузлами 57а, деј + К = і. 


2 7 
54" Sud Su“ 


Punc.: Дерево біноміального методу 
[] Seydel R.U. Tools for Computational Finance. 6th ed. / R.U. 
Seydel . - London: Springer-Verlag, 2017. - 508p. 


Біноміальний метод. Модель Кокса-Росса-Рубінштейна 


ə При у = 1 у дереві біноміального методу, кожні 2 кроки в 
часі повторюється початкове значення ціни на актив Sọ. 


0.45 т т т 
т Неее + + о + + 
0.4 ааа + + + 
к 
«НН ЕЛЕ + + + ЕЕ ЕЕ H 
зненненик юю + + + + жо ж 
0.353 анна а + + + + + + + + 
знане юю ++ + ++ + + 
+++++++ + + ++ + ++ 
зннененнею К+ + + + + + о +о+ ++ 
оз тнненачакта ж а + + + 
Ныне + я 
знанню + + + + + + +++ 
Неее + 
0.25 - НН + +++++ + Я 
Неее + + 
+++ жжжжяою 
+++ +++ + + + 
0.2} rrr + + 1 
+++ + + + +о+ + 
+++ + + 
+++ + 
0.15 +++++ 1 
+++ ++ + + + 
+++++ 
0.1 Чет 4 
s НН 
++++ 
+ 
0.05 ++++ 4 
++++ 
+++ 
++ 5 
0 | 1 1 1 
0 50 100 150 200 250 


Рис.: Дерево біноміального методу при М = 32 


[] Seydel R.U. Tools for Computational Finance. 6th ed. / R.U. 
Seydel . - London: Springer-Verlag, 22017. — 508p. 
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Біноміальний метод. Розв'язок рівнянь 


ә З рівнянь (4) та (10) отримуємо квадратне рівняння: 
и? 28 +7 =0, (13) 


де В = 1(уе—'^* + вк, 


ә Розв'язки: 


а= ти В – 402 - у, (14) 


и – а 
ә Розв'язки у випадку моделі Кокса-Росса-Рубінштейна: 


и = e7V^t + 0(/ 43), 


а= 1/и = 677V At, (15) 
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Біноміальний метод. Алгоритм 


ә Ініціалізація дерева. Прямий хід. 
ə 5 = So = 500 - корінь дерева 
ә Для і = 1,2,..., М обчислюємо: 


Sji = Soui di. (16) 


г Тепер вузли сітки (5; ;, &) - зафіксовані і в цих вузлах 

будемо обчислювати вартість опціону Vj; := У (5) 1, ti). 
ә Обчислення вартості опціону в вузлах. 

о В кінцевий момент часу tm вартість опціону У(5, їм) - 
відома. Тому вона відома і для кожного значення вартості 
акції 5; м = Sof dM, j = 0, ..., М. Це визначає значення 
У; м: Саї-опціон: 


М; м := тах(5; м - К,0). (17) 


Риї-опціон: 
М; м := тах(К - $; м, 0). (18) 
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Біноміальний метод. Алгоритм 


ә Обчислення вартості опціону в вузлах (продовження) 


ә Обернений хід. Обчилюємо значення У рекурентно для 
їм-1, їм-2, яз to: 


Му = е "Аруы + (1 - p) Vj ні). (19) 


ә Для європейського опціону значення Vo о обчислене за 
рекурентною схемою (19) і буде шуканим наближенням 
вартості опціону в даний момент часу 

ә У випадку американського опціону, опціон може бути 
виконаний в будь-який момент часу, тому схема (19) 
мусить бути модифікована таким чином: 

М пе = е "АУ + (1 = p) Мл) 
Мы = тах{ тах{ Sj; - К,0}, М2?" } для СаЙ-опщюну} 
Мы = тах{ тах К - 5,1,0}, Vpr} для Риї-опціону) 
(20) 
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Біноміальний метод. Алгоритм в цілому 


ә Обчислюємо At :- T/M. 


ә u,d,p - обчислюємо за формулами(14). 
500 =50 о 
5м = So oui dM, j =0,1,..., М 
Для американського опціону також обчислюємо 
5ы = Soo d'™ ana 0< i< M,j = 0,1,..., Ї. 

ә Обчислюємо М; м за формулами (17) - (18). 

ә Обчислюємо У; ; для і < М за формулами (19) для 
європейського опціону і (20) для американського. 

ә Отримане таким чином значення Voo - апроксимація ММ) 
шуканої вартості опціону М (50,0). 
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Метод скінченних різниць для моделі 


Блека-Шоулза 


Стельмащук Віталій Володимирович 


Львівський національний університет імені Івана Франка 
Кафедра інформаційних систем 


8 грудня 2022 р. 
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чиї 


Рівняння Блека-Шоулза 


ә Рівняння Блека-Шоулза: 


ду ыга дм 


Е 75 М =0 1 
Е + 252 1795 ' а) 
ә Рівняння Блека-Шоулза після заміни змінних (див. лекція 
б): 
T 2 
Yr = Уух, ХЕ (—оо,оо), тє (0, 75-) (2) 


ә Після заміни змінних, термінальні умови для ціни опціону 
М замінюються на початкові умови для у: 


у(х, 0) = (е(К+1)х/2 — е(к—1)х/2)+ при x Є (оо, оо); для Call 
у(х, 0) = (е(ко-1)х/2 — elkt1)x/2)+ при x Є ( -00, оо); для Put 
(3) 
к= (4) 
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Основи методу скінченних різниць. Різницеві 
апроксимації 


ә З розкладу в ряд Тейлора, кожна f Е С? задовольняє: 


по) = ПАННО ие) (5) 


де Є - деяка точка між x та х + Й. 


ө Рівномірна сітка: 


a ©. Хі < Хх Хіт S 
h:= X41 — Xi (6) 
fi := f (xi) 


ə Одностороння апроксимація: 


РО) = повно Оч) L olh vre С? (7) 
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Основи методу скінченних різниць. Різницеві 
апроксимації 


ə Двосторонні (центральні) апроксимації: 


Р(х) = Пон Ма O(h) Ее сЗ (8) 


(xi+1) — 28(хі) + Иж) 


Е 
Р(хі) = h2 


+ O(h?) yfe С" (9) 
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Основи методу скінченних різниць. Сітка 


® Двовимірна сітка вузлів для дискретизації задачі (2): 


Та? 
Дт := пра Вол Tmax === (10) 


Ах := (b — a)/m 


Тут ми замінюємо нескінченний інтервал —со < х < оо на 
деякий скінченний а = хуп < X < Xmax = Ё. 


е Значення а та Ь вибираються з умови: 
Ke? = Smin < $ < Smax = Ke? (11) 


ә Крім того, інтервал [а, b] має містити відповідні х-значення 
для So та К. 


@ w - розв'язок дискретизованої задачі, що відповідає задач! 


(2). 
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Основи методу скінченних різниць. Сітка 


Ах 


Рис.: Сітка 


П Seydel R.U. Tools for Computational Finance. 6th ed. / R.U. Seydel . - London: Springer-Verlag, 


2017. — 508p. 
е Позначення: 
Ty = „Ат для и = 0,1,..., Vmax 
хі := a + ГДх для і = 0,1,..., т 
Yiv := у(х}, Ty) 


АИ 
М/і и © Ури 


(12) 


ә уу; - визначені лише в вузлах сітки 
ә Похибка апроксимації е = Wi v — У; у залежить від вибору 
параметрів дискретизації итах, т, а, Б. 
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Явний (прямий) різницевий метод 


ә Апроксимації: 


дуг, Yiv+1 — У, 
, = ? 7 | А 
дт Ат OLAT) 


13) 
yiv Yi 1 = 2Yi v + Уї-1ш ( 
Й = Й р ? О А 2 
дх? Ах? PULAR 
ә Різницеве рівняння щодо и: 
Wiv+1 Муи Миа — 2Wiv + Мі у 
В y = 7 E РА 14 
Ат Ах? Се 
ә Позначення: Ж 
Е 
А = — 15 
Ах? ( ) 


ә Компактний запис: 


Wi +1 = Ан; ты + (1 — 2^) и + Аи и (16) 
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Явний (прямий) різницевий метод 


ә w”) := (miv, ...› Wmi) 
ə Матриця: 
1-22 А 0 2.5 0 
А 1-2^ à 5 0 
А := Ашы ПУ К Шш “е Re тб 
0 же А 1-2 ^ 
0 5 0 ^ 1-22 
(17) 
ө Різницеве рівняння в матричному записі: 
МОЇ) = Аш), у = 0,1,2,... (18) 
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Явний (прямий) різницевий метод 


т У+1 


Рис.: Шаблон явного різницевого методу 


[] Seydel R.U. Tools for Computational Finance. 6th ed. / R.U. Seydel . - London: Springer-Verlag, 


2017. — 508p. 


Явний різницевий метод не завжди стійкий. Умова стійкості: 


Ах? 
Ue Are (19) 
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Неявний (зворотній) різницевий метод 


ө Апроксимація похідної по часу: 


дуг ‚_ Yiv — Yiv- | 
s =e О(дт) (20) 


ә Компактний запис різницевої схеми: 
— Аи и + (1 + 2^) и = Аи зу = Мі у -1 (21) 
ə Матриця: 


1-2 -А 0 збу 0 
—À 1+2^ —^А не 0 
0 ы -А 1-22 -А 
0 еа 0 -А 1-22 
(22) 


А := Ањы У | 
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Неявний (зворотній) різницевий метод 


ә Різницеве рівняння в матричному записі: 
Аи!) = w”), у =0,1,2,....итах = 1 (23) 
ə На кожній ітерації методу необхідно розв'язувати СЛАР 
(23). Матриця - тридіагональна, тому складність 
розв'язування СЛАР - O(m). 


ә Схема неявного різницевого методу безумовно стійка 
УДт > 0. 


v+l 


i=l i іні 


Рис.: Шаблон неявного різницевого методу 


П Seydel R.U. Tools for Computational Finance. 6th eJ. / RU. Зеуаєї . - Гопаоп: Springer-Verlag, 
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Різницева схема Кранка-Ніколсона 


е Попередньо розглянуті різницеві схеми наближали ду з 


точністю О(Ат). Було б добре побудувати різницеву схему, 
яка є безумовно стійкою та наближає ду з точністю 
О(Ат?). 

ә Кранк та Ніколсон запропонували усереднити прямий та 
зворотній різницеві методи. В результаті отримуємо такий 
компактний запис різницевої схеми: 

СЕ Wiv _ 5 (мт — ЗМ u+ (24) 
+1, + Мі ЗМ уч + МіО уч) 


У+1 


Рис.: Шаблон різницевого методу Кранка-Ніколсона 
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Різницева схема Кранка-Ніколсона 


ә Матриці: 
1+А -8 0 0 
А | А 
5 1+А -3 0 
А := Acn := і zi (25) 
А А 
| o „| {ч | 
1-4 А 0 | 
А А 
à 1-Х А 0 
В:= Всем := |... о Ми (26) 
А А 
0 à 1-А à 
0 0 А 1-А 
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Різницева схема Кранка-Ніколсона 


ә Різницеве рівняння в матричному записі: 


Awt!) — Bw”), p= 0, 1, 2, ое Vmax — 1 (27) 


Теорема (Кранк та Нїколсон) 


ə Для у Є С" похибка наближеного розв'язку дорівнює 
О(Ат?) + О(Ах?) 

® Для кожної ітерації в часі м потрібно розв'язувати СЛАР з 
тридіагональною матрицею. 

ә Схема стійка для будь-якого Ат > 0 


Доведення: // Seydel R.U. Tools for Computational Finance. 6th 
ed. / R.U. Seydel . - London: Springer-Verlag, 2017. - 508p. 
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Алгоритм різницевої схеми Кранка-Ніколсона 


ә Вибираємо т, Vmax, Обчислюємо Ах, Дт. 


0 
® Обчислюємо wl = у(хі, 0) за даними початковими 


умовами (3), 0 < і < т. 


Формуємо матриці А та В. 


В циклі для и = 0,1,..., Vmax — 1 розв'язуємо СЛАР 


Awt!) = Bwl”) (28) 
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Крайові умови для європейських опціонів 


ө Крайові умови: 


у(х, т) = п(х, т)прих — —оо, 


у(х, т) — n(x, т)прих — оо. (29) 
ө Call опціон: 
n(x, T) = 0, 
ө Риї опцїон: 
= 1(к1)х+1(к—1)2т 
Е. | (31) 


ю(х,т) = 0. 
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Врахування крайових умов в схемі Кранка-Ніколсона 


ә Крайові умови в дискретизованій задачі виглядатимуть 
наступним чином: 


Ми Ии = п (а, ту), 


Ити = 2, Ту). (32) 
ә Вектор крайових умов для схеми Кранка-Ніколсона: 
п(а, тузі) + п(а, ту) 
дез ° (33) 
- 0 


2 (Ь, тил) + го, Ty) 


ә Матричний запис схеми Кранка-Ніколсона: 


Awt!) = Bw”) + d”), == 0, 1, 2, <.. Vmax 7 1 (34) 
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Моделювання відсоткових ставок. Процес 


Орнштейна-Уленбека. Модель Васічека 
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Стохастичний процес Орнштейна-Уленбека 


е Виділяють стохастичні процеси, в яких значення 
випадкової змінної в часі повертаються до її середнього 
значення (математичного сподівання). Такі процеси в 
англомовній літературі називають mean-reverting 
processes. 


® Процес Орнштейна-Уленбека: 


Тут 0 - коефіцієнт збігу випадкової змінної Х; до її 
середнього, и - середнє значення випадкової змінної Xt, 
с - волатильність (стандартне відхилення) випадкової 
змінної Ху, We - вінерівський процес. 
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Моделювання відсоткових ставок. Модель Васічека 


ə Насправді ринкова відсоткова ставка - це не стала 
величина, тому виникає потреба в моделюванні її 
поведінки. 

ə Модель Васічека передбачає, що ринкова відсоткова ставка 
г, слідує процесу Орнштейна-Уленбека з и > 0. 


ап = к(0 — гуаї + саму, (2) 


Тут к - коефіцієнт збігу випадкової змінної г; до її 
середнього, 9 > 0 - середнє значення випадкової змінної rt, 
с - волатильність (стандартне відхилення) випадкової 
змінної гу, И - вінерівський процес. 

ə Недоліком моделі є те, що вона не виключає від'ємних 
значень r. 
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Модель Васічека 


е Аналітичний розв'язок стохастичного диференціального 
рівняння моделі Васічека (виводиться з леми |то для е“): 


t 
r = пе “+0(1-е *") + гет | е“ ау; (3) 
0 


ә Змінна г, - нормально розподілена з математичним 
сподіванням: 


Е(һ) = пе + 0(1- е“) (4) 


та дисперсією: 
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Інші моделі відсоткової ставки 


ә Модель Кокса-Інгерсолла-Росса (Cox-Ingersoll-Ross model): 
айп = к(0 — гаї + oy redW: (6) 


виключає можливість негативних значень для змінної гу. 


ə Модель Халла-Вайта (Hull-White model): 


ап — (0+ = ок)дЕ + ordW: (7) 
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Облігації з нульовим купоном та оцінка їх вартості 


ə Облігація з нульовим купоном (zero-coupon bond) 
передбачає, що в певний момент часу Т у майбутньому 
інвестору виплатять її номінальну вартість (face value) X. 

ә Виникає питання: скільки грошей У варто заплатити 
інвестору за таку облігацію у кожен момент часу t < Т. 

о У випадку сталої відсоткової ставки г - це Present Value, 


тобто 
Y = X x eT), (8) 
ә У випадку ж змінної в часі відсоткової ставки /;: 
Y = X x eT "9%. (9) 


ә Значення r(s) визначаються стохастичними процесами 
моделей Васічека, Кокса-Інгерсолла-Росса чи Халла-Вайта. 
Замість аналітичного розв'язку для рівнянь з цих моделей 
можна використовувати метод Монте-Карло і 
усереднювати згенеровані траєкторії модельованого 
стохастичного процесу. 
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